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Low-Lying Eigenvalues of the Wilson-Dirac Operator * 
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An exploratory study of the low-lying eigenvalues of the Wilson-Dirac operator and their corresponding eigen- 
vectors is presented. Results for the eigenvalues from quenched and unquenched simulations are discussed. The 
eigenvectors are studied with respect to their localization properties in the quenched approximation for the cases 
of SU(2) and SU(3). 



1. Introduction 

In this contribution we consider the eigenvalues 
and corresponding eigenvectors of the hermitean 
operator 



Q = 75M/ (1 + 8k) , 



(1) 



where M is the Wilson-Dirac operator defined 
through 



+ {l+-Y M )Ul(x-n)1>(x-ii)} .(2) 

The hopping parameter k is related to the bare 
quark mass mo by k — (8 + 2m ) _1 . The gauge 
fields Uft(x) are mostly taken to be in SU(2) and 
for some cases in SU(3). The matrix Q is nor- 
malized in such a way that the eigenvalues A of 
Q satisfy < A 2 < 1. In the following we will 
concentrate on those eigenvalues of Q closest to 
zero, and refer to this part of the spectrum as the 
"low- lying" eigenvalues. 

The eigenvalues of Q are of crucial importance 
for understanding the performance of algorithms 
used in dynamical fermion simulations. In partic- 
ular, the scaling behaviour of the algorithms when 
approaching the chiral limit can be expressed in 
terms of the lowest eigenvalue of Q 2 Q| . 
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Furthermore, the eigenmodes themselves play 
an interesting physical role. In this work we will 
only address the space-time localization proper- 
ties of the eigenvectors belonging to low-lying 
eigenvalues. Additional more challenging ques- 
tions include their connection to chiral symme- 
try breaking and topology. Another important 
aspect is the contribution of the low-lying eigen- 
values to physical observables like the pion or p 
meson correlation functions. 

We compute the low-lying eigenvalues and the 
corresponding eigenvectors by using a modified 
Conjugate Gradient algorithm to minimize the 
Ritz functional 



(3) 



For the computation of higher eigenvalues the 
minimization of (||) is restricted to being or- 
thogonal to previously calculated eigenvectors. 
Convergence is accelerated by intermediate di- 
agonalizations in the space spanned by {^i} [||. 
This method has the virtue of being numerically 
stable and of providing rigorous error bounds 
for the eigenvalues calculated. Moreover the ap- 
proximate eigenvectors are obtained without ad- 
ditional cost. 

2. Low-lying eigenvalues for SU(2) 

We first investigate how the spectrum changes 
with the gauge coupling g 2 ,. Recall that for free 
Wilson fermions and k < 1/8, the lowest eigen- 
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Figure 1. Low- lying eigenvalues as a function of 
k for quenched SU(2) configurations at f3 = 2.3. 



Figure 2. Lowest eigenvalue on 8 4 (open sym- 
bols) and 16 4 (full symbols) lattices from un- 
quenched simulations as a function of K sca - 



value of Q 2 is 8-fold degenerate. Working in 
the quenched approximation on an 8 4 lattice at 
k = 0.12, we find that up to g$ = 4/ [3 around 0.8 
the lowest eight eigenvalues are nearly degenerate 
and well separated from the ninth lowest. Above 
g 2 ~ 0.8 the eigenvalues form a dense band. It 
is interesting to note that this value of f3, where 
the gap to the ninth eigenvalue closes, is roughly 
where one might expect the finite temperature 
phase transition. 

The k dependence of the n lowest eigenvalues 
(with n = 1,3,5,7,10) at /? = 2.3 is shown in 
Fig. 1 for a lattice size of 12 4 . The results in Fig. 1 
are averaged over 20 configurations each sepa- 
rated by 500 pure gauge updates. For k < 0.16 
the low-lying eigenvalues are almost degenerate 
and appear on top of each other in the graph. 
For k > 0.16 the eigenvalues drop significantly 
and their density becomes smaller. 

The coarse K-scan in Fig. 1 provides only an 
overall picture of the behaviour of the eigenval- 
ues. On individual configurations, used for Fig. 1, 
a much finer K-scan, guided by a perturbative ex- 
pansion in small changes in k, reveals almost-zero 



modes with values of A 2 that can drop to e.g. 
O(10~ 18 ), many orders of magnitude smaller than 
the lowest eigenvalue in Fig. 1. So far, it seems 
that these almost-zero modes appear above k c . If 
this observation is true for all coupling parameter 
values of interest, these almost-zero modes would 
be harmless. Otherwise they would indicate the 
breakdown of the quenched approximation since 
any occurance of such a mode could severely dis- 
tort the configuration sample. 

To investigate the situation for dynamical 
fcrmions we made unquenched runs for SU(2) at 
[3 = 2.3. In Fig. 2 we show the lowest eigen- 
value on 8 4 and 16 4 lattices. For the 16 4 lattice, 
the scan in the sea-K has a high resolution. On 
some of these configurations we also performed 
fine valence-K scans. In both cases we did not 
find any indication of the existence of the almost- 
zero modes. Fig. 2 also shows that for the 16 4 
lattice the lowest eigenvalue in unquenched simu- 
lations has a sharp minimum. This suggests that 
the lowest eigenvalue could be a good quantity to 
identify the phase transition. 
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Figure 3. Participation ratio for the eigenvalues 
in Fig. 1. 



Figure 4. 3-dimensional slices of p(x) at k = 
0.1475. The grey scale from bright to dark cor- 
responds to values of the density p(x) between 
and p max . 



3. Localization of the eigenmodes 

A simple quantity to study localization prop- 
erties of wave functions is the participation ratio 



where V is the lattice volume and (a, a) are Dirac .' - - 

and colour indices. For pointlike localization one . i 

ol (tains P — 0, whereas P — 1 for a uniform i 

distribution. The participation ratios which cor- : : : ... . ::: 

respond to the eigenvalues of Fig. 1 are shown in . 
Fig. ||. They dearly indicate strongly localized 

eigenmodes for n > 0.17, but do not give 1 a clear 

indication at lower k values. 

To investigate this question further we plot the 
local density p(x) — 1 1 \I/ (a^) 1 1 2 for 3-dimensional 
cuts through a lattice of size 16 3 48 for quenched 

SU(3) at 3 = 6.0. Fig. 4 indicates localization _,. _ _,. , , 

, £ _ .. . • - i figure 5. Same as bis. 4, but at k = 0.165. 

also lor k < n c « 0.157; however it is less pro- 
nounced than for n > k c (Fig. 5). 
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For a quantitative analysis of the localization 
we define 

Rv(c) = y#{ X '■ P{ X ) > C'Pmax} (5) 

which measures the relative number of lattice 
sites with a local density larger than a certain 
fraction c (0 < c < 1) of the maximum value 
Pmax- We find that practically all sites contribut- 
ing to Rv(c) are in a connected region, and that 
Rv(c) scales proportional to the inverse lattice 
volume. 

To characterize the localization we try for p{x) 
a parametrization p(x) ~ exp {— m; oc ||x||}, where 
|| a; || denotes the distance from the site where p = 
Pmax- The effective "localization mass" mi oc is 
the local rate with which p{x) decays. Taking for 
1 1 a; 1 1 the effective radius r g(c) — ^/2Rv{c)/tt 2 
one obtains m; oc = — ln(c)/r e ff. A plateau in the 
c-dependence of mi oc would then indicate that 
p(x) decays exponentially according to the above 
ansatz. 
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Figure 6. Localization mass on a 16 3 48 lattice 
from quenched simulations in SU(3). 



This is indeed observed as Fig. 6 shows for a 
value of k well below k c . One finds a nice sta- 
ble plateau in m/ oc completely consistent with an 
exponential localization. 



4. Conclusion 

We have presented an exploratory study of the 
low-lying eigenvalue spectrum of the hermitian 
Wilson-Dirac operator Q given in eq. ([!]). In the 
quenched approximation we find that the lowest 
eigenvalues drop sharply when the critical hop- 
ping parameter is reached. The isolated almost- 
zero modes we find in the quenched approxima- 
tion at [3 — 2.3 have not occurred in our un- 
qucnchcd runs. The eigenmodes corresponding 
to the low-lying eigenvalues show an exponential 
localization for k < k c . For k above k c they are 
even stronger (almost pointlike) localized. In gen- 
eral we did not see any qualitative difference be- 
tween the results from SU(2) and SU(3). 

Although the eigenvalues of Q can not be di- 
rectly related to the complex spectrum of M, 
we expect the qualitative behaviour of the eigen- 
modes to be similar. Of course, this relation has 
to be explored. It will also be interesting to relate 
the low-lying eigenvalues to topological properties 
of the theory and to quantify the contribution of 
the low-lying eigenvalues to physical observables. 
We hope to return to these questions in future 
publications. 
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